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Abstract 

From large basis set coupled cluster calculations and a minor empirical ad- 
justment, an anharmonic force field for silane has been derived that is con- 
sistently of spectroscopic quality (±1 cm~^ on vibrational fundamentals) for 



> 

cn 

o 

\^ . all isotopomers of silane studied. Inner-shell polarization functions have an 

o, 

O^ , appreciable effect on computed properties and even on anharmonic correc- 

Y^ . tions. From large basis set coupled cluster calculations and extrapolations 



to the infinite-basis set limit, we obtain TAEo=303.80ib0.18 kcal/mol, which 



c/3 

0^. includes an anharmonic zero-point energy (19.59 kcal/mol), inner-shell corre- 



lation (—0.36 kcal/mol), scalar relativistic corrections (—0.70 kcal/mol), and 
atomic spin-orbit corrections (—0.43 kcal/mol). In combination with the re- 
cently revised AiJ^o[Si(g)], we obtain Ai:f^_Q[SiH4(g)]=9.9±0.4 kcal/mol, in 
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between the two established experimental values. 



I. INTRODUCTION 



The spectroscopy and thermochemistry of the silane (SiH4) molecule have aroused in- 
terest from a number of perspectives. Its importance as a precursor for the chemical vapor 
deposition (CVD) of silicon layers has been discussed at length by Allen and Schaefer 0, 
who also review early theoretical work on the molecule. 

The spectroscopy of the tetrahedral group IV hydrides AH4 (A=C, Si, Ge, Sn, Pb) has 
been extensively studied. For a review of early work on AH4 (A=Si, Ge, Sn) the reader is 
referred to Ref. 0. 

A complete bibliography on experimental work on methane and its isotopomers would 
be beyond the scope of this work (see Refs. |Q,^ for detailed references): we do note that an 
accurate ab initio force field was computed by a team involving two of us. Based on this 
force field, a number of theoretical spectroscopic studies of the excited vibrational states of 
CH4 were recently studied: we note in particular a full-dimenstional variational study by 
Carter et al. @], a low-order perturbation theoretical/resonance polyad study by Venuti et 
al. 1^, and a high-order canonical Van Vleck perturbation theory study by Wang and Sibert 
0. We also note an accurate anharmonic force field on the isoelectronic NH^ molecule by 
two of us. 

The infrared spectrum of silane, SiH4, was first studied in 1935 by Steward and Nielsen 
P] and a set of fundamental frequencies for the most abundant isotopomer was first obtained 
in 1942 by Nielsen and coworkers. [^] 

The isotopomers of SiH4 have been the subject of considerable high-resolution experi- 
mental work; for instance, we note [0,0 for 28SiH4, 29SiH4, 30SiH4, [0-|T3 for ^sSiHgD, 



|T^,[T6| for ^'^SiHDa, and p!7| , p!8| for ^^SiD4. The molecule is of considerable astrophysical 
interest, having been detected spectroscopically in the atmospheres of Jupiter and Saturn 
| TP| and in the interstellar gas cloud surrounding the carbon star IRC-l-10 216 PU 



Until most recently, only fairly low-resolution data [^ were available for SiH2D2; as 
the present paper was being prepared for publication, a high- resolution study |^ of the 
{z/3, z/4, z/5, z/7, z/g} Coriolis resonance polyad appeared, in which assignments were facilitated 
by mixed basis set CCSD(T) and MP2 calculations of the quartic force field. 

One of the interesting features of the infrared spectra of silane is their pronounced local- 
mode character (e.g. [^), leading to complex resonance polyads. The strongly 'local' char- 
acter also inspired a study of the SiH4 spectrum up to seven quanta using algebraic methods 
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In the present work, we shall report a high-quality quartic force field that is of constant 
quahty for all the isotopomers of silane. A theoretical spectroscopy study by Wang and 
Sibert [^ is currently in progress on excited states and vibrational resonance polyads of 
SiH4 and isotopomers, using high-order (6th and 8th) canonical Van Vleck perturbation 



theory p6| and the force field reported in the present work. 

Since this can be done at very little additional computational expense, we shall also 
report a benchmark atomization energy and heat of formation of SiH4. The thermodynamic 
properties of silane are linked to a controversy concerning the heat of vaporization of silicon, 
which is of fundamental importance to computational chemists since it is required every 
time one attempts to directly compute the heat of formation of any silicon compound, be 
it ab initio or semiempirically. AifjQ[Si(g)] is given in the JANAF tables [^ as 106.6±1.9 



kcal/mol. Desai p8[ reviewed the available data and recommended the JANAF value, but 



with a reduced uncertainty of ±1.0 kcal/mol. Recently, Grev and Schaefer (GS) p9[ found 
that their ab initio calculation of the TAE of SiH4, despite basis set incompleteness, was 
actually larger than the value derived from the experimental heats of formation of Si((7), 
ii{g), and SiH4(5f). They concluded that the heat of vaporization of silicon should be revised 
upwards to Aifj_Q[Si(f7)]=108. 07(50) kcal/mol, a suggestion supported by Ochterski et al. 
||30|| . Very recently, however, Collins and Grev (CG) [^ considered the scalar relativistic 
contribution to the binding energy of silane using relativistic coupled cluster techniques 
within the Douglas-Kroll ||32| (no-pair) approximation, and found a contribution of -0.67 



kcal/mol. This would suggest a downward revision of the GS value of A_f/'jQ[Si(g)] to 107.4 
kcal/mol, which is in excellent agreement with a recent redetermination by Martin and Taylor 
p3| of 107.15±0.39 kcal/mol. (This latter value was derived by combining a benchmark ab 
initio calculation of the total atomization energy of tetrafluorosilane, TAEo[SiF4], with a 
very precise fluorine bomb calorimetric measurement |^| of Aiifj[SiF4(g)].) 



In addition, it was pointed out ^^ that the JANAF value of AifJ JSiH4(^)]=10.5±0.5 
kcal/mol is in fact the Gunn and Green [^ value of 9.5±0.5 kcal/mol increased by a correc- 
tion |3^ of +1 kcal/mol for the phase transition Si(amorphous)^Si(cr). (Gunn and Green 
considered this correction to be an artifact of the method of preparation and ignored it.) 

Clearly, a calibration calculation of TAEo[SiH4] might be desirable, and is the secondary 
purpose of the present study. Accurate thermochemical parameters of SiH4 (and other 
silicon compounds) are of practical importance for the thermodynamic and kinetic modeling 



of such processes as laser-induced chemical vapor deposition of silicon films from silane [^ 
the chemical vapor deposition of tungsten contacts for ULSI (ultralarge scale integrated 
circuit) chips by SiH4 reduction of WFg (e.g. PS|) and the generation of SiOxNy films by 



low-pressure chemical vapor deposition from mixtures of SiH4 with N2O and/or NH3 
(e.g. as antireflective coatings |^D| and for ultrathin capacitors ^^). (We also mention in 



passing the use of silane compounds in dentistry 

While GS's work was definitely state of the art in its time, the attainable accuracy for 
this type of compound may well have gone up an order of magnitude in the seven years 
since it was published: in a recent systematic study [^ of total atomization energies of a 



variety of first-and second-row molecules for which they are precisely known, procedures like 
the ones used in the present work achieved a mean absolute error of 0.23 kcal/mol, which 
dropped to 0.18 kcal/mol if only systems well described by a single reference determinant 
(as is the case with SiH4) were considered. In order to ascertain the utmost accuracy for 
hydrides, a zero-point energy including anharmonic corrections was found to be desirable 
P3| : this is obtained as a by-product of the accurate anharmonic force field which is the 
primary subject of the present contribution. 



II. COMPUTATIONAL METHODS 

All electronic structure calculations were carried out using MOLPRO 97 [Q running on 
DEC Alpha and SGI Origin computers at the Weizmann Institute of Science. 

The CCSD(T) [coupled cluster with all single and double substitutions (CCSD) ^5| sup- 
plemented with a quasiperturbative estimate of the contribution of connected triple excita- 
tions 1^] method, as implemented in MOLPRO |^, was used throughout for the electronic 
structure calculations on SiH4. For the Si('^P) atom, we employed the definition of Ref. | JS| 
for the open-shell CCSD(T) energy. 

The calculations including only valence correlation employed the standard Dunning cc- 
pVnZ (correlation consistent valence n-tuple zeta [^) basis sets on hydrogen and two dif- 
ferent variants of the cc-pVnZ or aug-cc-pVnZ (augmented cc-pVnZ [^,0) basis sets on 
Si. The first variant, cc-pVnZ+1, was used in the force field calculations, and includes an 



additional high-exponent d function ||52| to accommodate the greater part of the inner-shell 
polarization effect, which is known to be important for both energetic and geometric prop- 
erties of second-row molecules. ||5^ , |53[| The second variant, aug-cc-pV?2Z+2dlf |Q, includes 



two high-exponent d functions and a high-exponent / function, with exponents determined 
by successively multiplying the highest exponent already present for that angular momentum 
by a factor of 2.5. Such a set should give ||5^ an exhaustive account of the energetic effects 
of inner-shell polarization. 

Calculations including inner-shell correlation (not to be confused with inner-shell po- 
larization, which is an SCF-level effect) were carried out using the Martin- Taylor []54 



core 



correlation basis set. Relativistic effects were determined with the same basis set and as 
ACPF (averaged coupled pair functional |5^) expectation values of the first-order Darwin 



and mass- velocity operators [56, 



Optimizations were carried out by univariate polynomial interpolation. Force constants in 
symmetry coordinates were determined by recursive application of the central finite difference 
formula: the symmetry coordinates are defined in the same way as in previous studies [^0 



on the isovalent CH4 and NH4 molecules. The vibrational analyses were performed using a 



modified version of the SPECTRO program |58,59| running on an IBM RS6000 workstation 



at NASA Ames and the DEC Alpha at the Weizmann institute. The alignment conventions 



for the anharmonic constants of a spherical top follow the work of Hecht pO] and general 



formulae for these constants were taken from the paper by Hodgkinson et al. [BT|. Similar 



to previous work [§,^ on the spherical tops Be4 and CH4, the accuracy of the various 
spectroscopic constants was verified by applying opposite mass perturbations of ±0.00001 
a.m.u. to two of the hydrogen atoms, then repeating the analysis in the asymmetric top 
formalism. 

Finally, the reported zero-point energies include the Eq term [^ (which is the polyatomic 
equivalent of the a^ Dunham coefficient in diatomics) . 

III. RESULTS AND DISCUSSION 
A. Vibrational frequencies and anharmonic force field 

An overview of the basis set convergence of the computed bond distance, harmonic fre- 
quencies, and vibrational anharmonic corrections is given in Table 1. 

The effect of adding inner-shell polarization functions to the cc-pVTZ basis set is modest 
but significant (0.006 A) on the bond distance: the Si-H stretching frequencies, however, 
are affected by 20-25 cm^^. The bending frequencies are not seriously affected: somewhat 
surprising are the fairly strong effects on the vibrational anharmonicities (including, to a 
lesser extent, the bending frequencies). The overall behavior is in contrast to previous 



observations |£3[ for SO2 in which the inner-polarization effects on lower-order properties 
like geometry and harmonic frequencies are very noticeable but those on anharmonicities 
next to nonexistent, but is consistent with the very strong basis set sensitivity noted for the 



first three anharmonic corrections of the first-row diatomic hydrides by Martin p4 |. 



Likewise, a rather strong sensitivity with respect to basis set improvement from VDZ+1 
over VTZ+1 to VQZ+1 is seen for the Si-H stretching frequencies and all the anharmonic- 
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ities, even as the harmonic bending frequencies appear to be close to converged with the 
VTZ+1 basis set. It appears that in general, basis set sensitivity of anharmonicities of A-H 
stretches is much more pronounced than that of A-B stretches. 

The effect of inner-shell correlation, while nontrivial for the purpose of accurate calcu- 
lations, is quite a bit more modest than that of inner-shell polarization (as measured by 
comparing the cc-pVTZ and cc-p VTZ+1 results), and in fact is not dissimilar to what one 
would expect for a first-row molecule (e.g. CH4 [||). 

We will now consider computed fundamentals for the various isotopomers of silane with 
our best force field, CCSD(T)/cc-pVQZ+l. All relevant data are collected in Table 2. 

For ^^SiH4, ^^SiH4, and '^°SiH4, agreement between the computed and observed funda- 
mentals can only be described as excellent, with a mean absolute deviation of 2.5 cm~^. 
Agreement for the completely deuterated isotopomer ^^SiD4 is even better, with a mean 
absolute deviation of 1.9 cm~^. For the ^^SiHaD isotopomer, agreement is likewise excellent, 
with a mean absolute deviation of 2.1 cm^^. It would appear that the force field is certainly 
of good enough quality to permit assignments for the less well known isotopomers. 

For ^^SiHDs, the only precisely known bands are the Si-H stretch, z/i=2187.2070(10) 
1 ig, and the 1/5 degenerate bend, 850.680823(10) cm'^ fTg. Meal and Wilson [^, 



cm 



in their 1956 low-resolution study, assigned absorptions at 1573, 1598, and 683 cm"^ to 1/2, 
1^4, and z/g, respectively. Our calculations confirm this assignment and are on average within 
about 2 cm~^ of all the above bands, u^ was not observed by Meal and Wilson, and these 
authors speculated that it coincide with the 683 cm~^ (z/g) peak. Our own calculations pre- 
dict a splitting of about 5.7 cm~^ between z/3 and uq; B3LYP/VTZ+1 [^ infrared intensity 



calculations suggest that both bands should be observable. Inspection of the relevant spec- 
trum (Fig. 3 in Ref. [^) revealed that, at the resolution afforded by the equipment used, 
meaningful resolution between z/3 and uq becomes essentially impossible, especially given 
contamination (noted by Meal and Wilson) from a SiIl2D2 impurity with z/4=682.5 cm^^. 
Until most recently, the only available information for SiIl2D2 was the Meal and Wilson 



work. Our calculations, like those of Rotger et al. [|3l; unambiguously suggest assignment of 



the 1601 and 1587 cm~^ bands to V2 and v^, respectively, rather than the opposite assignment 
proposed by Meal and Wilson. We note that z/g is in a very close Fermi resonance with Z/5 + Z/9 
(the unperturbed levels being only about 10 cm^^ apart), despite the fairly small interaction 
constant /c569=-20.88 cm~^. Our calculations confirm the assignments for all other bands 
aside from vi and v^, which are calculated to be within 1 cm~^ of each other such that a 
meaningful decision on whether or not to exchange z/i and v^ is impossible. The Meal- Wilson 
empirical force field value of 844 cm^^ for v^ (which they were unable to observe) agrees well 
with our calculation as well as with the high-resolution value |2^ of 842.38121(9) cm~^. 



Of the very recent measurements by Rotger et al. p2[, all five bands in the Coriolis 
pentad (1/3, 1/4, z/5, 1/7, and i/g) are in excellent agreement with the present calculation (mean 
absolute deviation 1.1 cm~^). 

Among the sources of residual error in the quartic force field, neglect of inner-shell corre- 
lation and imperfections to CCSD(T) appear to be the potentially largest. As seen in Table 
1, inclusion of core correlation increases harmonic frequencies by as much as 7 cm~^ in this 
case. The effect of correlation beyond CCSD(T) was seen to work in the opposite direction 
for the first-row diatomic hydrides |]B^; in the present work, we have compared FCl/VDZ-l-1 
and CCSD(T)/VDZ-|-1 harmonic frequencies for the SiH diatomic in the X ^11 and a ^E~ 
states, and found a reduction in Ue of 4 and 10 cm~^, respectively. (The FCI-CCSD(T) 



difference for u,, was found in Ref. 0] to converge very rapidly with the basis set.) Since 
FCI frequency calculations in a reasonable-sized basis set for SiH4 are simply not a realistic 
option, we have taken another track. 

We have assumed that the computed CCSD(T)/VQZ-|-1 force field is fundamentally 
sound, and that any residual error would mostly affect the equilibrium bond distance and the 
diagonal quadratic force constants. We have then taken our quartic force field in symmetry 
coordinates, substituted the computed CCSD(T)/MTcore bond distance (which agrees to 
four decimal places with the best experimental value), and have iteratively refined the four 
diagonal quadratic force constants such that the four experimental fundamentals of ^^SiH4 
are exactly reproduced by our calculation. The final adjusted force field is given in Table 3 



and is available in machine-readable format from the corresponding author. 

As seen in Table 2, our computed fundamentals for the other isotopomers with the ad- 
justed force field are in essentially perfect agreement with experiment where accurate values 
are available. Discrepancies arise for some modes of SiH2D2, SiHDs, and SiD4 where only 
low-resolution data are available. Particularly the discrepancy for V2 of SiD4 is completely 



out of character: the experimental difficulties involved in its determination fl^ suggest that 
perhaps the experimental value may be in error. (A discrepancy of 1.2 cm~^ for V2 in SiHsD 
is halved upon accounting for a Fermi resonance 2v^ ~ z/2-) We hope that our computed 
force field will stimulate further spectroscopic work on SiH4 and may serve as a basis for 
studies employing more sophisticated vibrational treatments, such as the variational tech- 
niques variational techniques very recently applied to methane |^ or high-order canonical 
Van Vleck perturbation theory. As noted in the Introduction, a study of the latter type is 
already in progress. p5| 

B. Geometry 

At the CCSD(T)/MTcore level, we compute a bond distance of 1.4734 A, which we 
know from experience []66| , |67| should be very close to the true value. Ohno, Matsuura, 
Endo, and Hirota (OMEHl) |6^ estimate an experimental r^, bond distance of 1.4741 A 
without supplying an error bar; in a subsequent study (0MEH2) ||69|, the same authors, using 
two different methods, obtain 1.4734(10) A ("method I") and 1.4707(6) A ("method 11"), 
respectively, where uncertainties in parentheses are three standard deviations. The deviation 
between the (diatomic approximation) "method 11" value and our present calculation is more 
than an order of magnitude greater than usual for this level of ab initio theory, while the 
"method I" value agrees to four decimal places with our calculation. (Normally, because 
of neglect of correlation effects beyond CCSD(T) which have the tendency [^ to lengthen 
bonds by 0.0002-0.0006 A, we expect our computed bond distance to be slightly short, rather 
than too long.) The computed bond distance of Rotger et al., 1.4735 A at the CCSD(T)[all 



electron] level in a mixed basis set which does not contain any core correlation functions, is 
likewise in excellent agreement with the 0MEH2 "method I" value. 

C. Atomization energy of SiH4 

Using a 3-point geometric extrapolation A+B.C~^ from the SCF/AVr2Z+2dlf (r2=T,Q,5) 
atomization energies, we find an SCF limit component of the total atomization energy of 
259.83 kcal/mol, only marginally different from the directly computed SCF/AV5Z+2dlf 
value of 259.82 kcal/mol and only 0.05 kcal/mol larger than the GS result. 

The CCSD valence correlation component was extrapolated using the 2-point formula 
|7| A + B/n^ from AVnZ+2dlf (ra=Q,5) results; thus we obtain a CCSD limit of 64.26 
kcal/mol, which is 0.8 kcal/mol larger than the largest basis set value (63.45 kcal/mol) 
and 1.4 kcal/mol larger than the largest basis set value of GS (62.86 kcal/mol). Using the 
alternative 3-point extrapolation [^ A + B/{1 + 1/2)*^ from AVnZ+2dlf (r2=T,Q,5) we 



obtain a somewhat smaller basis set limit of 63.92 kcal/mol; however, as discussed in Ref. 
this procedure appears to systematically underestimate basis set limits and was found 



7^ to yield excellent agreement with experiment largely due to an error compensation with 
neglect of scalar relativistic effects. 

At 0.81 kcal/mol, the extrapolated basis set limit contribution of connected triple exci- 
tations is quite modest, and differs by only 0.02 kcal/mol from the largest basis set value 
of 0.79 kcal/mol. In fact, it is largely immaterial whether the extrapolation is done from 
AVr7,Z+2dlf (ra=T,Q) or from AVr2Z+2dlf (ra=Q,5), and we obtain essentially the same re- 



sult for the (T) contribution as GS (0.82 kcal/mol). This is an illustration of the fact ^3 
that connected triple excitations generally converge more rapidly with basis set than the 
CCSD correlation energy. 

Adding up the two basis set limit values, we find a valence correlation component to TAE 
of 65.05 kcal/mol; given the essentially purely single-reference character of the SiH4 wave 
function there is little doubt that the CCSD(T) limit is very close to the full CI limit as well. 
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As noted by GS, the contribution of inner-shell correlation of SiH4 is negative: we find 
-0.365 kcal/mol compared to their -0.31 kcal/mol. The spin-orbit contribution is trivially ob- 
tained from the Si(^P) atomic fine structure [0 as -0.43 kcal/mol, while our computed scalar 
relativistic contribution, -0.70 kcal/mol, is essentially identical to the CG value. Finally, we 
obtain TAEe=323.39 kcal/mol. 

The anharmonic zero-point vibrational energy (ZPVE) from our best force field (includ- 
ing Eq) is 19.59 kcal/mol. This is very close to the value of 19.69 kcal/mol obtained by GS 
as an average of estimated fundamentals and CISD/TZ2P harmonic frequencies: the com- 
putational effort involved in improving this estimate by a mere 0.1 kcal/mol would therefore 
have been hard to justify if the anharmonic force field would not have been required for an- 



other purpose. Also, from past experience |75], we know that such good agreement between 



rigorous anharmonic ZPVEs and estimates cannot be taken for granted for hydrides. 

Our best TAEe and ZPVE finally lead to TAEo=303.80 kcal/mol, to which we attach an 
error bar of about 0.18 kcal/mol based on previous experience [^. This should be compared 
with the GS largest basis set result of 303.03 kcal/mol (or 302.36 kcal/mol after applying the 
CG scalar relativistic contributions) or the value derived from JANAF heats of formation of 
Si(g), H(g), and SiH4(g), 302.62 kcal/mol. 

If we consider alternative values for Aif;o[Si(g)] of 108.1±0.5 kcal/mol (GS), 107.4±0.5 
kcal/mol (applying CG to the latter value), or 107.15±0.38 kcal/mol (Martin & Taylor 
||33|), we would obtain from our calculation AHjQ[Si}l4{g)] values of 10.8±0.5, 10.1±0.5, and 
9.9±0.4 kcal/mol, respectively. Only the first of these values cannot be reconciled with Gunn 
and Green; the very similar values derived from the CoUins-Grev-Schaefer and Martin- Taylor 
Ai7jQ[Si(g)] agree to within accumulated error bars with both the JANAF and Gunn-Green 
values for the heat of formation of silane. While our best value of 9.9±0.4 kcal/mol at 
first sight slightly favors the Gunn-Green value (in which the Si{cr)—>-Si{amorph) transition 
enthalpy [^] was considered an artifact of the manner of preparation), the difference is "too 
close to call". We contend that our calculated value is more reliable than either experiment. 
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IV. CONCLUSIONS 

From accurate ab initio calculations and a minor empirical adjustment, a quartic force 
field for silane has been derived that is consistently of spectroscopic quality (±1 cm^^ on 
vibrational fundamentals) for all isotopomers of silane studied here (^^SiH4, ^^SiH4, '^°SiH4, 
^^SiHsD, ^^SiH2D2, ^^SiHDs, and ^^SiD4). As in previous studies on second-row molecules, 
we found that inner-shell polarization functions have an appreciable effect on computed 
properties, and for hydrides this apparently includes the vibrational anharmonicities. 

From large basis set coupled cluster calculations and extrapolations to the infinite-basis 
set limit, we obtain TAEo=303.80±0.18 kcal/mol, which includes an anharmonic zero- 
point energy (19.59 kcal/mol), inner-shell correlation (—0.36 kcal/mol), scalar relativis- 
tic corrections (—0.70 kcal/mol), and atomic spin-orbit corrections (—0.43 kcal/mol). In 
combination with the recently revised Aif?Q[Si(g)], 107.15±0.39 kcal/mol [Q, we obtain 
AifjQ[SiH4(g)]=9.9±0.4 kcal/mol, intermediate between the JANAF and Gunn-Green val- 
ues of 10.5±0.5 and 9.5±0.5 kcal/mol, respectively. 
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TABLES 
TABLE L Basis set convergence of computed bond distance (A), harmonic frequencies (cm^^^ 
and anharmonic corrections (cm^^) of ^^SiH4; effect on inner-shell correlation. 







VDZ 


VDZ+1 


VTZ 


VTZ+1 


VQZ+1 


MTcore 


MTnocore 


Te 




L49076 


1.48572 


1.48504 


1.47952 


1.47872 


1.47339 


1.47736 


UJl 




2242.0 


2249.3 


2225.7 


2250.3 


2262.7 


2270.6 


2264.1 


UJ2 




978.3 


982.4 


983.1 


985.5 


983.4 


991.4 


987.1 


W3 




2253.1 


2259.6 


2227.9 


2254.7 


2266.5 


2275.4 


2268.3 


UJ4 




925.8 


933.6 


932.5 


933.8 


930.8 


937.2 


935.3 


^1 




2167.0 


2175.6 


2154.8 


2174.3 


2185.0 






V2 




965.6 


970.0 


964.9 


969.0 


968.3 






vz 




2173.1 


2181.0 


2155.3 


2175.1 


2185.2 






Z/4 




912.2 


920.3 


913.6 


917.9 


915.1 






Wi - 


- v\ 


74.96 


73.65 


70.90 


75.96 


77.73 






ijJ2 - 


- V2 


12.77 


12.43 


18.14 


16.51 


15.10 






•^3 - 


-Vy, 


80.04 


78.54 


72.58 


79.60 


81.25 






LOa, - 


- ^A 


13.65 


13.34 


18.96 


15.89 


15.70 







The CCSD(T) electron correlation method has been used throughout. 
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TABLE II. Comparison of computed and observed fundamentals (cm ^) for isotopomers of 
silane. 



i 


Vi 


L^i 


i^i 


Wi -Vi 


Vi 




CCSD(T)/ 


best 


best 


best 


Expt. 




cc-pVQZ+1 


adjusted 


adjusted 


adjusted 










28SiH4 






1 


2185.0 


2264.2 


2186.9 


77.34 


2186.873254(80) |k ] 


2 


968.3 


986.0 


970.9 


15.10 


970.93451(6) [0 


3 


2185.2 


2270.1 


2189.2 


80.84 


2189.189680(66) |k ] 


4 


915.1 


929.1 


913.5 

29SiH4 


15.68 


913.46871(4) [[![ 


1 


2184.9 


2264.2 


2186.8 


77.39 


2186.8281(5) [0 


2 


968.3 


986.0 


970.9 


15.09 


970.94856(22) |^ | 


3 


2183.7 


2268.4 


2187.7 


80.70 


2187.6494(1) [0 


4 


913.8 


927.8 


912.2 

30SiH4 


15.63 


912.18278(8) [0 


1 


2184.9 


2264.2 


2186.8 


77.43 


2186.7855(6) [0 


2 


968.3 


986.0 


971.0 


15.08 


970.95790(110) ||li| 


3 


2182.2 


2266.8 


2186.2 


80.57 


2186.1963(1) [0 


4 


912.6 


926.6 


911.0 

28SiH3D 


15.59 


910.97921(12) ||i| 


1 


2184.9 


2265.7 


2187.4 


78.38 


2187.40066(5) Q 


2 


1590.7 


1630.9 


1592.8(a) 


38.08(a) 


1593.9595(10) Q 


3 


914.5 


928.5 


912.9 


15.60 


912.991(1) [[4| 


4 


2184.4 


2270.0 


2188.4 


81.51 


2188.50418(4) Q 


5 


949.2 


966.1 


950.6 


15.48 


950.576(1) [[4 


6 


784.6 


795.4 


784.2 

28SiH2D2 


11.29 


784.324(1) [[4 


1 


2184.3 


2269.9 


2187.7 


82.19 


2189 [0 


2 


1579.4 


1621.1 


1581.3 


39.83 


1587 [0 (b) 


3 


942.0 


958.4 


942.7 


15.69 


942.74106(4) Q 


4 


681.3 


689.6 


681.3 


8.29 


681.62394(3) || 


5 


840.1 


854.0 


842.3 


11.71 


842.38121(9) [§|| 


6 


1597.2 


1640.6 


1599.7 


40.87 


1601 


7 


861.1 


873.4 


859.5 


13.84 


859.750104(4) Q 


8 


2183.6 


2267.2 


2187.3 


79.81 


2183 |2l) ?? 


9 


743.7 


753.0 


742.4 

28SiHD3 


10.65 


742.64029(3) [g 


1 


2183.1 


2268.5 


2186.6 


81.90 


2187.2070(10) ^ 


2 


1570.8 


1611.4 


1572.4 


39.01 


1573 


3 


676.4 


683.9 


675.2 


8.72 


[682] O (c) 
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4 


1596.0 


1640.5 


1598.7 




41.77 


1598 


5 


850.2 


863.4 


850.6 




12.84 


850.680823(10) ||I6| 


6 


682.1 


690.6 


682.4 

26 


SiD4 


8.14 


683 [|i| 


1 


1562.6 


1601.7 


1563.8 




37.84 


1563.2(10) [L7| 


2 


687.9 


697.5 


689.8 




7.75 


685.2(2) |i7 


3 


1595.2 


1640.5 


1598.0 




42.48 


1598.44919(43) (I||, 1598.45(5) 


4 


675.3 


682.6 


674.1 




8.48 


674.2(15) 



(a) If Fermi resonance 2z/8 ^ V2 is accounted for (2z/|=1563.9 cm~^, A;288=21.393 cm"^, 
1/2 = 1587.1 cm^^) we obtain z/2=1594.6 cm"^, and 2z/8=1556.5 cm"^. 

(b) in fact doublet at 1584 and 1591 cm"^; we suggest assignment of 1584 cm~^ to z/2 and 
of 1591 cm~^ to possibly v^ + z/g 

(c) not observed; valence force field estimate. Authors of Ref. [^ speculate that it coincides 
with the 683 cm~^ band. 
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TABLE III. Quadratic, 


cubic and quartic 


force constants (aJ/ 


'A'^radian") 


for SiH4 


Fu 


3.04428 


-^22 


0.41777 


F44 


2.92753 


F7A 


-0.08914 


F77 


0.51105 


i^iii 


-6.72559 


F22I 


-0.16483 


F441 


-6.55397 


i^741 


0.05353 


F77I 


-0.14909 


F222 


-0.02902 


-^662 


-0.11954 


Fg62 


0.13392 


-F992 


-0.29002 


-^654 


-6.51514 


F954 


-0.07640 


-f984 


-0.03587 


-^987 


0.41927 


Fun 


12.58584 


-?^2211 


-0.00037 


-^4411 


12.80551 


Fyau 


0.05928 


-P77II 


-0.06795 


-^^2221 


-0.01703 


Fm2i 


0.04882 


-^9621 


-0.00005 


-^9921 


0.09065 


-^6541 


12.78806 


-f9541 


0.02617 


^9841 


0.05749 


-^9871 


-0.11162 


-^2222 


0.16525 


-^6622 


-0.10547 


-^6633 


-0.22839 


-f9622 


-0.08371 


-^9633 


-0.02617 


-^9922 


0.17090 


-f9933 


0.30175 


-F9542 


0.18309 


-^8762 


-0.04530 


F4444 


13.10285 


-^5544 


12.93973 


F7444 


0.03663 


-f8544 


0.05928 


F7744 


0.00877 


-^8754 


-0.00357 


-f8844 


-0.21095 


F7774 


-0.24594 


-^8874 


-0.08029 


-^7777 


0.26859 


-f8877 


0.75864 
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